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Evidence  for  surface  a particle  clusters  in  ““‘Ag  and  *^Au 

from  the  (e,a)  reaction 


W.  R.  Dodge''’ 

8200  Raymond  Lane 
Potomac,  Maryland  20854-3727 


The  “‘Ag  and  ‘^Au  (e,p)  and  (e,a)  energy  and  angular 
distributions  were  measured  at  6 electron  bombarding  energies 
between  50  and  115  MeV  at  the  National  Institute  of  Standards  and 
Technology.  The  “‘Ag  and  *^Au  (e,p)  angular  distributions  exhibit 
an  asymmetric  component  which  increases  from  0 to  50%  as  the  proton 
energy  increases  from  the  proton  Coulomb  barrier  to  26  MeV.  The 
“‘Ag  and  ‘^^Au  (e,a)  angular  distributions  exhibit  an  asymmetric 
component  which  increases  from  0 to  30%  as  the  a energy  increases 
from  the  a particle  Coulomb  barrier  to  26  MeV.  We  conclude  that 
the  asymmetric  component  of  the  proton  and  a particle  yields  are 
the  result  of  direct  or  semidirect  processes  rather  than  resonance 
processes  and  hence,  because  of  the  short  mean  free  paths  of  a 
particles  in  nuclear  matter,  give  evidence  for  the  existence  of  a 
particle  clusters  in  the  nuclear  surface. 


^The  experimental  work  reported  in  this  paper  was  carried  out  while 
the  author  was  an  employee  of  the  National  Institute  of  Standards 
and  Technology. 


/ 


' ,v  '.,■  ■ 'i::^S 


tv ' 


• • ’,!“  ’ ' - . ' \ 'Si,-'  .''...I  (,  ■ ' ’•  ' . '.''  V''..-  ■ ii*''' 

"'-U''.'''=^^  .,'';J;''f  fa-,, ’4i/’'>  «fSp^Sli 

;',.  i-.f  -f-  =:  ,?  k ^ ^ -S,  --'^W  ^ ^ 

‘SV'"''-'  f v:.;"  o-f 

' - 


J^t'-  -t':  .:■>^  ;* 


Vr*  ^ W 


■i/<m  s - 4-  ,>:r '%-, '^®|#,Ss,fft 


t 'V.Tt,t.  ^ ' '**?^*f >t';  -.W^'  f 

■■  ^",1  ,'i-,’---’^', 


■ '.. ' ,•  ".■R.'.j  . ',  ’ ' ’i.  '. 

■-'  '-iM''  ' ''  ' ' 


' ''"'  ,' -V; t4'' '‘i 


? Vi- ; 


) ■* 


Ai: 


..  y. 


I . INTRODUCTION 


The  hypothesis  that  nucleons  in  the  surface  of  a heavy  nucleus 
are  correlated  in  a particle  clusters  of  the  type  hypothesized  by 
Wilkinson’  has  received  experimental  support  from  several  different 
sources.  These  sources  are  the  ratio  of  the  theoretical  and 
experimental  values  of  the  reduced  widths  for  a decay  in  a 
emitters,’  knockout  reactions  produced  by  protons  and  alphas,  (a, 
2a)  coincidence  experiments,^  and  K'  meson  capture  from  circular 
orbits.’  To  these  traditionally  quoted  sources  of  evidence  for 
surface  a particles  has  recently  been  added  evidence  from  a 
transfer  reactions  such  as  (d,^Li)  (n,a)  scattering  experiments,'* 
(p,p'a),^  and  fast  tt  nuclear  scattering.^  These  experiments  have 
all  contributed  to  the  plausibility  of  the  hypothesis  of  nuclear 
surface  a particles.  Of  these  experiments,  the  (a, 2a)  and  (p,p'a) 
coincidence  experiments  are  the  most  convincing,  while  the  fast 
N(77',7r')N*  inelastic  scattering  experiments,  in  which  the  product 
nuclei  were  detected  by  their  radioactivity,  can  be  explained  in 
terms  of  theories  which  do  not  assume  a clusters  according  to  Chang 
et  al.^ 

Arguments  for  the  occurrence  of  a particles  in  the  nuclear 
surface  are  based,  in  part,  on  the  fact  that  in  the  nuclear  surface 
the  nucleon  density  is  low,  and  hence  the  Pauli  principle  does  not 
inhibit  formation  of  nucleons  into  a clusters  as  strongly  as  it 
does  in  the  nuclear  interior.  Brink  and  Castro  have  studied  the 
importance  of  a clustering  in  the  interior  regions  of  heavy  nuclei 
by  studying  the  existence  of  such  phenomena  in  infinite  nuclear 
matter.®  They  showed  that  at  normal  nucleon  densities  a uniform 
nucleon  density  distribution  is  more  stable  than  the  a cluster 
nucleon  density  distribution,  indicating  that  a types  of 
correlations  are  not  important  in  the  interior  regions  of  heavy 
nuclei.  However,  at  nucleon  densities  of  a. 3 of  the  central 
nuclear  density  of  .068  nucleons/ fm^,  a transition  occurs  and  the 
a cluster  distributions  becomes  more  stable  than  the  uniform 
density  distribution,  indicating  that  a clusters  should  be  very 
important  in  the  surface  of  a heavy  nucleus.  While  these  a 
particle  clusters  are  of  obvious  importance  in  a description  of  the 
texture  of  the  nuclear  surface  and  as  examples  of  nuclear 
correlations,  their  relevance  to  theories  of  photonuclear  physics 
was  not  clear.  Recently,  however,  Kuo,  Blomqvist,  and  Brown’,  KKB, 
have  pointed  out  that  the  presence  of  a clusters  in  the  nuclear 
surface  could  partially  resolve  a discrepancy  between  Brown's 
microscopic  model  of  the  photonuclear  effect  and  experiment  in  ^°®Pb. 
According  to  KKB,’  these  alphas  would  tend  to  limit  the  volume  in 
which  protons  oscillate  against  neutrons  to  the  nuclear  interior, 
since  protons  and  neutrons  are  bound  together  strongly  in  the  a 
particle.  Furthermore,  these  alphas  could  provide  a basis  for  the 
assumption  of  the  macroscopic  hydrodynamic  models  that  the  proton 
and  neutron  relative  velocity  is  zero  at  some  nuclear  radius.’ 
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Motivated  by  these  theoretical  conjectures  about  the  existence 
of  a clusters  in  the  nuclear  surface,  we  have  searched  for  effects 
in  the  (e,a)  reaction  which  could  be  attributed  to  these  clusters. 
We  measured  the  (e,p)  and  (e,a)  energy  and  angular  distributions 
from  Be,  “‘Ag  and  ’^^Au  at  several  electron  bombarding  energies 
between  50  and  115  MeV.  Defining  the  asymmetry  of  the  measured 
angular  distribution  to  be  (f-b)/(f+b),  where  f and  b are  equal  to 
da/dn  integrated  from  0 to  7r/2  and  7r/2  to  tt,  respectively,  our  “‘Ag 
and  ’’^Au  da(e,p)/df2  exhibit  an  asymmetry  which  increases  from  a 0 
to  50%  as  the  proton  energy  increases  from  the  Coulomb  barrier  to 
25  MeV.  In  contrast,  our  da(e,a)/dn  have  an  asymmetry  which 
increases  from  ~ 0 to  3 0%  as  the  a energy  increases  from  the 
Coulomb  barrier  to  25  MeV.  Below  the  Coulomb  barrier,  our 
da(e,a)/dn  are  backward  peaked.  Reactions  which  exhibit  angular 
distribution  asymmetries  as  a function  of  the  particle  kinetic 
energy  are  typical  of  direct  reactions. Therefore,  we  conclude 
that  the  asymmetry  of  our  da(e,a)/dn  must  be  attributed  to  a direct 
or  semidirect  process  involving  alphas  which  are  not  formed  by  slow 
statistical  processes  during  de-excitation  of  the  nucleus. 


II.  EXPERIMENTAL  PROCEDURE 

The  National  Institute  of  Standards  and  Technology's  (formerly 
the  National  Bureau  of  Standards)  linear  electron  accelerator, 
linac,  provided  a spatially  resolved  and  momentum  analyzed  beam  of 
electrons  which  was  incident  on  self-supporting  “‘Ag,  *^Au,  and  Be 
targets.  The  “‘Ag,  ’^Au,  and  Be  target  mass  thicknesses  were  3.10, 
5.23,  and  14.95  mg/cm^,  respectively.  A magnetic  spectrometer"  was 
used  to  provide  particle  momentum/ charge  analysis  and  five  Si 
surface  barrier  detectors  of  thicknesses  varying  from  150  to  400 
and  positioned  in  the  spectrometer  focal  plane  provided  either 
particle  energy  or  energy  loss  which  was  used  in  particle 
identification.  In  addition,  three  of  the  focal  plane  detectors, 
FPDs,  were  backed  up  by  500  /xm  surface  barrier  detectors.  Above  ~ 
3.25  MeV,  protons  were  transmitted  through  the  FPDs.  Hence,  the 
proton  energy  loss  peak  in  the  FPDs  was  displaced  from  the  total 
absorption  position  to  a position  below  the  triton  peak,  t,  for 
most  of  the  proton  energies  measured  in  this  experiment.  Between 
« 3.25  and  8.2  MeV  protons  and  alphas  could  be  simultaneously 
detected  in  both  the  focal  plane  and  backing  detectors.  At 
energies  above  8.2  MeV,  protons  and  alphas  could  not  be 
simultaneously  detected  in  the  focal  plane  detectors  because  of  the 
finite  pulse-height  acceptance  of  the  data  logging  analogue-to- 
digital  converters,  ADCs,  and  hence  above  8.2  MeV,  protons  were 
only  counted  in  the  backing  detectors.  The  relative  detection 
efficiencies  of  the  focal  plane  and  backing  detectors  above  8.2  MeV 
were  compared  and  found  to  be  .96T.01,  .95^.01  and  .97T.01. 

The  momentum  acceptance  interval,  Ap,  of  each  detector  was 
determined  with  ^‘°Po  and  *'‘®Gd  a sources  (a  energies  of  5.304  and 


2 


3.18  MeV) . The  momentum  acceptance  interval  was  corrected  for  the 
fact  that  the  focal  plane  detectors  were  round  and  the  electron 
beam  spot  width  projected  on  the  target  was  not  equal  to  the  a 
source  width.  This  correction,  ®(W,,Wb) , to  Ap  was 


(1) 


where  W,  and  Wg  are  the  width  of  the  a source  and  electron  beam, 
respectively,  and  and  are  the  ratios  of  these  quantities  to 
the  focal  plane  detector  diameter  times  the  magnet  magnification. 
This  correction  to  Ap  was  never  larger  than  6%.  Data  from  each 
were  normalized  relative  to  the  central  FPD  by  comparing  the 
energy  spectra  of  each  detector  in  energy  intervals  .92  MeV  wide 
where  the  spectra  overlapped. 

A typical  pulse  height  spectrum  produced  by  the  charged 
disintegration  products  of  Be  is  shown  in  Fig.  1.  The  energy  of 
the  alphas  and  protons  is  4 . 5 MeV.  The  continuous  background  shown 
between  the  prominent  peaks  produced  by  the  charged  particle 
disintegration  products  shown  in  Fig.  1 had  two  components:  (1)  a 
contribution  from  momentum  analyzed  particles  produced  in  the 
target  ladder,  (2)  a continuous  contribution  from  electron,  photon, 
and  neutron  induced  events  not  momentum  analyzed  and  not  produced 
in  the  target  or  target  ladder.  Background  from  (1)  was  measured 
with  the  target  ladder  retracted  and  was  found  to  be  negligible. 
Information  about  background  (2)  was  provided  in  regions  of  the 
pulse  height  spectra  in  regions  between  well  resolved  mass  groups. 
Background  component  (2)  was  only  important  for  the  highest  energy 
protons  where  the  energy  deposited  in  the  rear  detector  became  < 1 
MeV.  In  these  cases,  background  (2)  was  subtracted  from  the  pulse 
height  spectra  by  a non-linear  least  squares  fitting  program.  The 
counts  extracted  from  the  various  mass  groups  in  the  FPD  pulse 
height  distributions  after  (2)  were  subtracted  were  divided  by  the 
number  of  incident  electrons,  e,  by  the  mean  energy  of  the  energy 
acceptance  interval,  2Ap/p,  of  the  detector  appropriate  to  the 
particle  charge  state  and  mass  to  correct  for  the  spectrometer 
magnet  dispersion. 

The  data  were  corrected  for  energy  loss  in  the  target  using 
energy  loss  data  compiled  by  Janni.*^  The  mean  target  thickness 
assumed  in  the  energy  loss  calculations  was  half  the  actual  target 
thickness  divided  by  the  cosine  of  the  angle  between  the  plane  of 
the  target  foil  and  the  spectrometer  magnet  symmetry  axis.  The  a 
energy  loss  in  the  “‘Ag,  ’’^Au,  and  Be  targets  never  exceeded  .15, 
.16,  and  1.5  MeV,  respectively. 

The  other  components  of  the  apparatus  used  in  the  experiment 
(i.e.,  the  linac  and  electron  beam  momentum  analyzing  and  handling 
system)  have  been  described  previously."  The  incident  electron 
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beam  momentum  transmission  interval,  (Ap/p)^,  was  ±1%.  The  incident 
electron  beam  current  was  measured  with  two  ferrite  toroid 
induction  current  monitors.  Beam  current  measurements  made  with 
these  monitors  deviated  from  each  other  by  < 1.5%  during  the 
experiment.  The  ferrite  monitors  were  calibrated  against  a 
standard  water  cooled  NIST  beam  dump.  The  calibration  constant  did 
not  depend  on  energy  between  20  and  110  MeV.  Typically,  all  data 
at  a given  angle  were  taken  in  a 12  hour  period. 

Energy  spectra  were  obtained  for  “‘Ag,  ‘’“^Au,  and  Be  as  a 
function  of  the  angle,  between  the  incident  electron  beam  and 
the  detected  particle,  and  as  a function  of  the  bombarding  electron 
energy  Eq.  Data  for  “‘Ag  and  *^Au  were  taken  at  = 34°,  48.,  62°, 

90°,  and  146°  with  Eq  = 50  and  110  MeV.  In  addition,  data  for  “‘Ag 
and  ‘^Au  at  34°  were  obtained  for  Eq  = 80  and  110  MeV,  and  at  90°  for 
Eq  = 75  MeV  and  100  MeV.  Data  for  Be  at  34°  were  obtained  at  = 
50,  80,  100,  and  110  MeV  and  at  90°  at  50,  75,  100  and  110  MeV. 


III.  RESULTS 

The  proton  and  alpha  yields/MeV/target  nuclei/sr/e  at  the 
various  angles  for  a given  value  of  Eq  were  least  squares  fitted  to 
an  expansion  in  Legendre  polynomials.  The  coefficients  of  the 
angular  distribution  expressed  in  the  center  of  mass  system,  cms, 
are  given  by 


A,  = + rc(c  + l) 


'1  + 1 


^1-1 


(2C  + 3)  (2<!-l) 


with  r = 


2T^ 


c 1. 


(2) 


In  this  expression  q is  the  virtual  photon  momentum  transfer,  Mp  and 
Mj  are  the  particle  and  target  mass  numbers,  and  Tp  is  the  particle 
kinetic  energy.  The  momentum  transfer,  q,  is  not  measured  in  this 
experiment,  but  for  Ag  and  Au  the  maximum  possible  momentum 
transfer,  2Eo/c,  results  in  a difference  between  most  angular 
distribution  coefficients  expressed  in  the  laboratory  or  cms  within 
errors  of  the  coefficients.  The  coefficient  Aq  is  shown  in  Fig.  2 
and  the  ratios  of  coefficients  A,/Ao,  Aj/Aq,  and  Aj/A^  are  shown  in 
Figs.  3 through  13  for  “‘Ag  and  ‘’^Au.  The  curve  for  Aq  multiplied! 
by  471  represents  the  energy  spectrum  integrated  over  particle 
emission  angles.  The  uncertainties  in  the  input  data  have 

contributions  from  counting  statistics  and  from  the  effective 
target  thickness.  The  uncertainties  in  the  effective  target 

thickness  are  generated  by  electron  beam  heating  changes  in  the 
angle  between  the  incident  beam  and  the  target  foil.  In  cases 
where  the  angle  between  the  normal  to  the  plane  of  the  target  foil 
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and  the  incident  electron  beam  was  small,  the  uncertainty  in  the 
effective  target  thickness  was  the  dominant  uncertainty.  The  ratio 
of  the  coefficients  Aj/A^  and  A3/A0  were  least  squares  fitted  to  an 
expression  of  the  form  ao+ajTp  over  certain  energy  intervals  with 
statistically  significant  x*  The  angular  distribution  asymmetries 
(defined  in  the  INTRODUCTION)  derived  from  these  fits  together  with 
the  average  asymmetries  computed  in  the  given  range  of  Tp  and 
computed  in  an  energy  interval  10  MeV  wide,  starting  where  the 
asymmetries  become  positive,  are  given  in  Table  I.  Fig.  14  shows 
the  differential  energy  spectrum  of  the  Be(e,a)  reaction  at  34°  for 
Eq  = 50  and  110  MeV. 

The  salient  features  of  the  energy  and  angular  distributions 
are:  (1)  the  “‘Ag  and  *’^Au(e,a)  angular  distribution  asymmetries 
increase  by  « 30%  as  Eq  is  increased  from  50  to  110  MeV,  (2)  the 
“‘Ag  and  *^Au(e,p)  angular  distribution  asymmetries  increase  by  ~ 
50%  as  Eq  is  increased  from  50  to  110  MeV,  and  (3)  the  ’’^Au(e,a) 
total  energy  distribution  peak  moves  to  a slightly  lower  energy  as 
Eq  increases  from  50  to  110  MeV.  In  contrast,  the  “‘Ag(e,a)  and 
(e,p)  and  the  ‘^Au(e,p)  energy  distribution  peaks  remain  fixed  in 
energy  as  Eq  changes  from  50  to  110  MeV. 

The  “‘Ag  and  '^Au  proton  and  alpha  yields/MeV/target 
nuclei/ sr/e  as  a function  of  the  incident  electron  energy,  Eq  at  34° 
and  90°  are  shown  in  Figs.  15  through  18.  The  data  were  averaged 
over  somewhat  larger  particle  energy  intervals  than  were  data 
presented  in  the  angular  distributions  in  the  interest  of 
compactness.  The  excitation  data  shown  in  Figs.  15  through  18  can 
be  converted  into  differential  cross  sections  by  the  virtual  photon 
difference  unfolding  method. 


^ do  (k,  6) 
dCl 


N(E^  ,k)  W(£„  ,k) 

f — - / -^dk 

three  three 


(3) 


where  [Y(Eoi)]e  is  the  yield/MeV/target  nuclei/target  nuclei/sr/e  at 
an  angle  of  6 and  N(Eo,)  is  given  by 


N{E^,k) 


ft 

2 

f 

1 + 

sin 

71 

k •^0  j 

1 

2E,(E,-k) 

m^k 


2 (E^-k)  1 
^0  J 


(4) 


Cross  sections  for  *’^Au(e,p)  and  (e,a)  derived  with  these 
assumptions  are  shown  in  Fig.  19.  Although  the  uncertainties  are 
large  and  the  points  few,  a significant  difference  is  seen  between 
the  trend  of  the  (e,p)  and  the  (e,a)  da/dn.  The  same  remarks  apply 
to  the  “‘Ag(e,p)  and  (e,a)  da/dn  shown  in  Fig.  20. 


5 


IV.  DISCUSSION 


The  primary  motivation  for  this  experiment  was  to  search  for 
anomalies  in  the  yields  of  a particles  from  Ag  and  Au  at  electron 
bombarding  energies  above  50  MeV  which  might  be  related  to  a 
clusters  in  the  nuclear  surface.  Since  Be  is  known  to  exhibit 
considerable  a clustering,  the  Be(e,a)  yields  were  intended  to 
provide  a reference  frame  with  which  to  compare  our  “‘Ag  and 
‘^Au(e,a)  yields.  By  anomalies  we  mean  features  of  the  (e,a)  energy 
and  angular  distributions  which  can  not  be  explained  in  terms  of 
pre-equilibrium  (precompound)  or  equilibrium  statistical  models. 
The  (e,Q:)  yields  could  be  due  to:  (1)  quasi-free  electron 
scattering  from  literal  a clusters  with  the  ejection  of  recoil  a 
cluster  from  the  nucleus;  (2)  absorption  of  a virtual  photon  by  the 
usual  photonuclear  photon  absorption  mechanisms,  including  photon 
absorption  by  a quasideuteron,  followed  by  precompound  and 
eventually  compound  nucleus  formation  and  evaporation  of  an  alpha; 
(3)  photon  absorption  by  a quasideuteron  with  the  fast  (p,n)  pair 
interacting  with  the  remaining  nucleons  via  the  same  or  similar 
interactions  responsible  for  (p,a)  and  (n,a)  scattering;  (4)  direct 
interaction  of  a virtual  photon  with  an  a cluster  and  the  remainder 
of  the  nucleus  in  the  sense  of  the  direct  photodisintegration 
models.  This  list  is  not  exclusive  but  includes  mechanisms  we 
believe  most  likely  to  be  responsible  for  the  observed  a yields. 
Items  2,  3 and  4 involve  conventional  photonuclear  processes. 

Assuming  the  existence  of  a clusters,  it  is  probable  that  all  of 
the  above  items  would  contribute  to  our  (e,a)  yields.  If  a certain 
fraction  of  our  a yields  result  from  direct  processes,  such  as 
outlined  in  items  1,  3,  and  4,  the  existence  of  "surface"  a 

clusters,  because  of  the  short  mean  free  paths  of  alphas  in  the 
nuclear  interior,  is  established.  On  the  other  hand,  the  compound 
nuclear  or  statistical  models  of  item  2 have  no  predictive  powers 
concerning  the  origin  or  existence  of  a clusters  since  all 
dynamical  aspects  of  the  evaporation  process  are  implicitly 
contained  in  the  reaction  or  inverse  cross  section  data  which  are 
an  input  into  the  calculation. 

We  now  consider  the  possibility  that  a significant  fraction  of 
our  a yields  are  produced  by  the  mechanism  of  item  1,  namely, 
quasifree  scattering.  Because  of  the  kinematic  incompleteness  of 
an  (e,a)  experiment  (i.e.,  the  scattered  electron  and  the  a were 
not  detected  in  coincidence)  , these  data  can  only  indicate  the 
plausibility  of  such  a process,  but  cannot  unambiguously  indicate 
its  existence.  However,  our  (e,a)  and  (e,p)  measurements,  together 
with  an  a particle  energy  loss  in  the  target  comparable  to  the 
proton  energy  loss  in  an  (e,e'p)  experiment,  indicate  that  the 
(e,e'a)  counting  rate  in  a heavy  nucleus  should  be  » .05%  of  the 
(e,e'p)  counting  rate.  Hence,  (e,e'a)  coincidence  experiments  are 
still  far  in  the  future  even  at  laboratories  with  electron  linacs 
with  100%  duty  factors.  Furthermore,  few,  if  any  (e,a)  experiments 
have  been  reported. 
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Here  we  estimate  the  magnitude  of  the  quasifree  a scattering 
component.  If  this  component  were  produced  by  elastic  scattering 
on  '‘He,  then  the  angular  distribution  of  these  alphas  in  the 
laboratory  system  would  be  given  by  the  Mott  scattering  cross 
section 


da 

dQ 


^ Eq  ) cos^d^  1 +2p  + p^sin^0^ 


(5) 


where  is  the  angle  between  the  elastically  scattered  a and  the 
incident  electron  of  total  energy  Eq,  p = Eq/M^,  M„  is  the  a particle 
mass,  ^(q^)  is  the  elastic  scattering  form  factor  of  '‘He.  If  a real 
nucleus  has  "a  clusters",  this  angular  distribution  will  be 
modified  by  the  momentum  distribution  of  these  alphas,  by 
scattering  of  these  alphas  by  other  nucleons  in  the  nuclear 
surface,  by  the  Coulomb  barrier,  and  by  the  number  of  "a  particles" 
with  mean-free  paths  such  that  these  alphas  escape  from  the 
nucleus . 


An  estimate  of  the  total  cross  section  for  recoil  a-production 
from  a '‘He  target  with  a energies  above  a cutoff  energy,  T^,  can 
be  obtained  by  integration  of  the  Mott  cross  section,  neglecting 
final  state  interactions  (assuming  only  the  Coulomb  interaction  is 
important)  and  assuming 

^(g2)  = g-<r2><7V3^  (6) 


which  is  a good  approximation  for  q^  < Here  q 

is  the  usual  3 -vector  momentum  transfer.  The  total  elastic 

electron  scattering  cross  section  on  '‘He, 


^^0  j I I 3 


2<r2> 


In 


(1  ^^2p)x- 


(1  +p) 


-—EQ^<r^> ^ 

3 (l  + p)2-p2^2 


1 

1 ll 

{l  + p)^-p^x^ 

[(l-Hp)2-p2x2]2jJ 

(7) 


where 


= cos^dj 


max 


T ■ (1  +p)  ^ 

^min  ^ ^ ' 


It  ■ + 

^ ' ^min 


2M„) 


(8) 
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The  total  number  of  recoil  alphas  from  elastic  electron 
scattering  on  '‘He  as  a function  of  Eq  is  shown  in  Fig.  21  for 
various  values  of  The  number  of  a particles  to  be  associated 
with  quasifree  scattering  depends  on  the  minimum  kinetic  energy 
which  an  a must  have  for  the  probability  for  penetration  through 
the  Coulomb  barrier  to  be  « 1.  Brink  and  Castro's*  infinite  nuclear 
matter  calculation  indicated  that  the  Fermi  momentum  at  nucleon 
densities  where  a clustering  becomes  important  is  about  200  MeV/c 
(kp  ~ 1 fm‘)  corresponding  to  a kinetic  energies  of  5 MeV.  On  the 
other  hand,  Dubost  ^ assumed  an  average  kinetic  energy  of  10 
MeV  for  a particles  in  the  nuclear  surface  in  analyzing  their  (p,a) 
data.  Igo,^  following  Rasmussen, on  the  basis  of  an  optical  model 
analysis  of  (a,a')  data,  conjectured  that  a particles  might  have  a 
stable  existence  where  the  nuclear  density  is  a 5%  of  the  central 
nucleon  density.  For  ‘^Au  this  corresponds  to  a radius  of  8 fm  and 
a Coulomb  barrier  energy,  T^,  of  14  MeV.  In  order  to  obtain  an 
estimate  of  the  number  of  quasifree  (e,a)  events  we  would  observe 
if  N„  a particles  existed  in  a real  nucleus,  we  use  the  momentum 
distribution  of  the  oscillator  well**^ 


Pa  (P) 


1 


4 

3 


/ 

V 


(9) 


The  oscillator  well  parameter  p^  = 198 . 7 (Z/A) MeV.  The  momentum 
distribution  of  a nucleon  in  the  oscillator  well  after  scattering 
is  given  by  6^„(P)  , where  P^=p^+q^-2p*q.  This  momentum  distribution 
is  to  be  integrated  over  all  orientations  and  magnitudes  of  p 
compatible  with  p > Pc  = Vm^Tc  . The  number  of  quasifreely 
scattered  particles,  N„,  with  energies  above  the  Coulomb  barrier  Tc 
is 

N,  = I (ff")  «.(gr»)  d(cose.)  , dO) 


where 

Jf  P,(P)p^cfpdQpg 

Pig 

is  the  fraction  of  recoil  particles  with  momenta  P >Pc.  The 
expression  9ta(q)  for  the  oscillator  well  is  derived  in  Appendix  I. 

While  a Fermi  gas  model  is  not  appropriate  for  an  a particle 
gas,  we  have  also  computed  9?F(q)  for  a Fermi  gas  momentum 
distribution  in  the  limit  of  infinite  dilution. 


8 


In  this  case, 
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for  kp  + qs^Pc-^q, 
St^iq)  = 


(12) 


' q'' 

\ ^FJ 
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16 
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for  kp  = Pc>q- 


Note  that  the  expression  for  Rp  = > q is  the  response  function 
given  by  deForest  and  Walecka’^  for  quasifree  (e,p)  scattering.  The 
Fermi  momentum  kp  was  taken  to  be  315  (Z/A)*'^  MeV.  For  *^Au  with  6^^ 
= 66.1°,  corresponding  to  a minimum  a recoil  energy  of  1 MeV  at  Eo 
= 110  MeV,  Eq.  (10)  predicts  .8  x 10'^°  N„  for  an  oscillator  momentum 
distribution  with  p^  = 147  MeV/c,  p^  = 376  MeV/c,  and  .4  x 10'^°  N„ 
for  a Fermi  momentum  distribution  with  = 232  MeV/c.  For  “‘Ag, 
the  corresponding  estimates  are  .30  x 10’"  N„  and  .25  x 10'^^  with 
Pc  = 318  MeV/c.  Increasing  p„  and  kp  by  a factor  of  1.5  increases 
the  ^’^Au  estimates  to  . 6 x 10'^^  N„  and  .9  x 10'^^  N„,  while  changing 
Tjnin  to  3 MeV  decreases  the  Au  estimate  by  a factor  of  10.5.  Thus 
estimates  obtained  from  Eq.  (10)  are  subject  to  large  uncertainties 
because  of  their  sensitivity  to  changes  in  T^  and  p„,  and  also 
because  the  impulse  approximation  may  not  be  valid  at  these 
electron  energies.  Of  course,  only  the  asymmetric  component  of  our 
“‘Ag(e,a)  and  *’^Au(e,a)  energy  spectra,  integrated  over  a energies 
and  angles,  can  be  attributed  to  quasifree  scattering.  At  Eq  = 110 
MeV,  these  are  13  and  14%,  respectively,  of  the  total  yields. 
Taking  into  account  theoretical  uncertainties,  and  using  p^  = 185 
MeV/c  and  T^  = 2 MeV,  our  data  for  Eq  = 110  MeV  are  compatible  with 
2.3±2i  and  .5±^4  alphas  in  the  nuclear  surface  of  Ag  and  Au, 
respectively.  Correction  for  final  state  interactions  would 
increase  these  numbers  significantly. 

Recoil  protons  from  quasifree  electron  scattering  should  also 
contribute  to  our  “‘Ag  and  ‘’^Au  110  MeV  proton  energy  spectra. 
Figure  22  shows  the  total  number  of  protons  elastically  scattered 
above  a cutoff  energy  T^  as  a function  of  Eg,  as  predicted  by  Eq. 
(10)  appropriately  modified.  Since  the  proton  recoil  energy  is 
four  times  as  large  as  the  a recoil  energy,  and  the  Coulomb  barrier 
is  half  as  large,  recoil  protons  should  be  clearly  discernible. 
Information  on  the  proton  momentum  distribution  from  (e,e')  and 
(e,e'p)  already  exists  for  *^C.  Inelastic  (e,e')  data‘^  show  the 
width  of  the  '^C  quasifree  elastic  scattering  peak  is  a 32  MeV  for 
Eq  = 98  MeV  and  6^  = 135°.  Hence,  the  width  of  the  quasifree 
recoil  proton  energy  distribution  for  Eq  = 110  MeV  will  have  a 
contribution  from  the  Fermi  momentum  distribution  of  qAq/Mp  ~ 6 MeV 
and  will  be  centered  at  Ep  = 18  MeV  for  0p  = 34°.  However,  unlike 
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alphas  in  which  are  unbound,  the  *^Au  proton  separation  energy 
is  5.8  MeV.  The  total  number  of  recoil  protons  above  the  Coulomb 
barrier,  which  is  « 9.4  MeV,  plus  the  proton  separation  energy  is 
.7  X 10'^*  protons/Au/e  using  % of  Eg.  (11)  with  Pc  = 15.2  MeV.  The 
*^Au(e,p)  data  shown  in  Fig.  5 have  only  a slight  inflection  at  18 
MeV  which  contains  a small  fraction,  <20%,  of  the  total  proton 
yield  of  1.4  x 10'^^  protons/Au/e.  Moreover,  this  inflection  is  also 
seen  in  the  50  MeV  *^^Au  data  at  the  same  proton  energy  and  is  not 
seen  in  the  ‘“‘Ag(e,p)  data  shown  in  Fig.  11.  Furthermore,  both  our 
“‘Ag  and  *^^Au  (e,p)  yields  fit  virtual  photon  isochromats  derived 
using  the  assumption  that  only  the  proton,  scattered  electron,  and 
a nucleus  in  the  ground  state  are  in  the  final  state  to  within  20%. 
Hence,  while  proton  shell  binding  energies  and  other  data  needed 
for  a proper  impulse  approximation  calculation  are  not  well  known 
for  electron  energies  of  110  MeV,  corresponding  to  reduced  electron 
wavelengths  of  1.8  fm,  and  while  the  impulse  approximation  is 
certainly  better  for  (e,p)  than  for  (e,a)  reactions,  our  failure  to 
observe  a measurable  quasifree  peak  in  our  proton  yields  leads  us 
to  discount  quasifree  scattering  as  a viable  mechanism  with  which 
to  explain  our  (e,a)  yields. 


We  now  consider  the  characteristics  of  the  (e,a)  energy  and 
angular  distributions  which  would  be  produced  by  the  mechanisms 
suggested  in  item  2 listed  at  the  beginning  of  this  section.  The 
photonuclear  version'®  of  the  statistical  model  predicts  that  the  a 
yields,  Y„,  are  given  by 


^0 


thres 


rjk)  NiE^.k) 

'£r,(k)  k 


(13) 


The  sum  over  the  partial  widths  FjCk)  is  over  all  open  decay 
channels.  Here,  (k)  is  the  total  photon  absorption  cross  section 
at  a photon  energy  k,  N(Eo,k)  is  defined  by  Eq.  (4),  and 

r^(Jc)  « J T p{k-S^-A-T)  dT 

thres 

In  Eq.  (14)  ^^^(T)  is  the  reaction  cross  section  for  the  ith 

particle  of  mass  number  Aj  on  a nucleus  of  mass  number  (A  - Aj)  , and 
p is  the  density  of  states  in  the  daughter  or  residual  nucleus.  In 
the  Fermi  gas  model, 

p(k-S,-A-T)  = eV*'*-'.--*-"  (15) 

where  a is  the  level  density  parameter  and  is  generally  taken  to  be 
a function  of  the  charge  and  mass  of  the  residual  nucleus,'’  S;  is 
the  separation  energy  for  the  ith  particle,  and  T,n«x  defined  by  j 
Eo”Sj-A-T,^  = 0.  The  pairing  energy.  A,  is  22 (Vh.  for  an  even-odd  and 
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44/\/a  for  an  even-even  nucleus. 


The  cross  section/MeV  , da„/dT, 


„/dT„,  is  given  by 


(16) 


Figs.  23  and  24  show  the  da/dT„  obtained  from  Eg.  (16)  for  ‘^Au 
(e,p)  and  (e,a)  reactions  using  the  values  of  a,  A,  and 
suggested  by  Dostrovsky  ^ ai-*’  The  photon  absorption  cross 
section  was  approximated  by  a Lorentz  function  which  fit  the 
^^Au(7,n)  data  in  the  energy  regime  of  the  giant  dipole  resonance 
and  by  the  quasideuteron  model  of  Levinger^^  with  the  Levinger 
parameter , if,  equal  to  10.6.  The  quenching  factor,  which  determines 
how  fast  the  giant  resonance  cross  section  is  turned  off  and  how 
fast  the  quasideuteron  cross  section  is  turned  on,  was 
Decreasing  the  level  density  parameter,  a,  from  14.9/MeV  to  10/MeV 
increases  the  cross  section  by  a 20.  Obviously  a value  of  the 
level  density  parameter  between  these  values  would  fit  the  '^Au(e,p) 
data.  However,  the  same  set  of  parameters  grossly  overestimates 
the  *^Au(e,p)  cross  section.  The  poor  quality  of  the  fits  reflects 
the  general  inadequacy  of  this  simple  version  of  the  statistical 
model.  Good  fits  to  (7,0:)  energy  spectra  obtained  at  lower  photon 
endpoint  energies^^  are  the  result  of  renormalization  of  the  cross 
section  magnitudes  or  use  of  unrealistic  level  density  parameters, 
a,  coupled  with  the  cutoff  of  the  bremsstrahlung  photon  spectrum. 
Within  the  framework  of  this  model,  the  d^/dn„/dT„  are  isotropic, 
although  Ericson  and  Strutinski^^  modified  the  model  so  that  angular 
distributions  symmetric  about  90®  can  be  obtained.  In  general, 
this  model*®  predicts  a differential  cross  section  d^a  (e, a) /dn„dT„ 
which  reaches  a maximum  a few  MeV  above  threshold  of  the  reaction 
cross  section  Coulomb  barrier  and  then  decreases  rapidly. 

Gabriel  and  Alsmiller,^  GA,  have  modified  the  conventional 
statistical  model  by  incorporating  an  intranuclear  cascade  model 
and  photon  absorption  by  the  quasideuteron  mechanism^*  to  provide 
both  a mechanism  for  increasing  the  cross  section  at  photon 
energies  above  the  giant  resonance  and  for  introducing  a forward 
asymmetry  in  the  emitted  neutron  and  proton  angular  distributions. 
A still  more  sophisticated  version  of  GA's  calculation  which 
includes  photon  absorption  by  tt  photoproduction  on  single  nucleons 
has  been  used  by  Adler  ^ al^  to  interpret  the  yield  of  alphas 
produced  by  500  MeV  bremsstrahlung  on  *^Au.  The  GA  calculation, 
modified  by  Adler  et  al^  to  include  the  effects  of  Coulomb  barrier 
penetration,  predicts  an  isotropic  a angular  distribution  which  is 
entirely  isotropic  in  the  center  of  mass  system.  This  isotropic 
angular  distribution  transforms  to  the  laboratory  system  to  da/dn,^ 
« aoPo(cos^)  + a,Pi(cos^),  where  aj/ao  « .054  ~ However,  the 

angular  distributions  of  Adler  ^ al^  exhibited  a forward  peaked 
anisotropic  component  of  16%  of  the  total  yield  which  lead  them  to 
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conclude  this  fraction  of  their  yield  resulted  from  direct 
processes . 

The  da/dT/unit  energy  loss  by  radiation,  0^^^,  where 


(17) 


is  referred  to  as  the  differential  cross  section/equivalent  quantum 
by  Adler  ^ al^^.  If  we  convert  our  data  to  a similar  quantity 
using  for  N(Eo,k)  the  virtual  photon  spectrum  defined  in  Eg.  (4), 
we  obtain  for  Au  a .51  mb/eq.  quantum  at  Eq  = 110  MeV  at  the  peak 
of  the  energy  spectrum  which  can  be  compared  to  2.2  mb/eq.  quantum 
obtained  by  Adler  ^ al^^.  The  differential  cross  sections/eq. 
quantum  obtained  by  Adler  ^ al^  and  by  this  work  are  shown  in  Fig. 
25.  The  experimental  and  theoretical  curves  of  Adler  ^ al^^  have 
been  renormalized  to  agree  with  our  peak  cross  section.  According 
to  the  latter's  calculation  of  the  energy  distribution  of  photons 
leading  to  the  emission  of  at  least  one  a,  our  a yield/eq.  quantum 
at  Eq  = 110  MeV  should  be  about  8%  of  their  yield/eq.  quantum  at  500 
MeV.  A linear  extrapolation  of  our  excitation  data,  shown  in  Figs. 
14  through  17  to  Eq  = 500  MeV  should  be  ~ 18%.  The  differences 
between  our  linear  extrapolation  and  the  calculation  of  Adler  et 
al^^  can  be  understood  in  terms  of  the  energy  distribution  of 
photons  leading  to  the  emission  of  at  least  one  a which  has  a 
dramatic  increase  (due  to  the  A resonance)  above  the  photopion 
threshold. 

Several  pre-equilibrium  or  precompound  statistical  models^ 
have  been  developed  which  aspire  to  provide  a unified  description 
of  the  statistical  properties  of  an  excited  nucleus  at  all  times 
after  an  exciting  event,  including  time  intervals  short  compared  to 
those  envisaged  in  the  statistical  models’*’^'*’^  discussed  above.  In 
these  models  the  states  of  the  system  are  described  in  terms  of  the 
number  of  particles  and  holes  they  contain  and  state  densities  are 
usually  evaluated  assuming  equally  spaced  single  particle  states. 
Complex  particles  are  fabricated  by  two-body  single  nucleon  pickup 
processes.  Early  versions  of  this  model  by  Cline, which  treated 
complex  particle  emission  in  the  framework  of  Griffin  and  Blann,^ 
ignored  the  Pauli  principle,  shell  structure,  and  angular  momentum. 
Later  calculations  by  BSt^k  ^ el^^  can  predict  the  relative 
integrated  energy  spectra  of  p,  d,  t,  and  a in  a scattering 
experiments^®  to  well  within  an  order  of  magnitude,  thus 
establishing  its  superiority  over  evaporation  statistical  models 
for  describing  the  gross  features  of  proton  and  a particle  induced 
complex  particle  energy  spectra. 

Neglect  of  angular  momentum  by  pre-equilibrium  models  has 
restricted  the  predictive  powers  of  this  model  to  the  isotropic 
component  of  complex  particle  yields.  We  speculate  here  on 
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difficulties  the  pre-equilibrium  models  may  have  after  angular 
momentum  is  considered.  Whatever  mechanism  is  responsible  for  our 
(S/P)  yields,  at  least  three  more  two-body  single  nucleon  pickup 
processes  are  necessary  to  produce  an  (e,a)  event.  While  single 
nucleon  pickup  processes  are  forward  peaked,  in  general,  after 
three  consecutive  such  interactions  needed  to  fabricate  an  a 
particle,  the  emitted  a particle  angular  distribution  would  be  less 
forward  peaked  than  the  angular  distribution  of  the  particle  which 
initiated  the  cascade.  However,  our  (e,p)  and  (e,a)  angular 
distributions  are  both  forward  peaked  and  not  significantly 
dissimilar.  The  same  remarks  apply  to  (q:,p)  and  (a, a)  angular 
distributions  at  an  incident  a energy  of  58  MeV^°.  These  facts 
imply  that  the  hypothesis  of  sequential  two-body  interactions  will 
have  to  be  augmented  to  satisfactorily  account  for  the  experimental 
data.  Of  course,  even  with  the  inclusion  of  angular  momentum, 
asymmetric  angular  distributions  will  not  be  predicted  by 
statistical  models  because  of  the  basic  statistical  model 
assumption  of  random  phase  relations  between  matrix  elements. 

Moreover,  the  asymmetry  does  not  seem  to  have  a trivial 
kinematic  origin,  since  transformation  of  an  isotropic  angular 
distribution  to  the  laboratory  system  fails  by  an  order  of 
magnitude  to  provide  the  asymmetry  observed  in  the  laboratory  for 
both  the  (a,p) , (a, a)  and  our  (e,p) , (e,a)  experiments.  Indeed,  if 
we  assume  that  our  angular  distributions  can  be  explained  by  a pre- 
equilibrium statistical  model  with  a few  of  the  nucleons  in  the 
nucleus  participating  in  the  reaction  (a  typical  situation  shortly 
after  the  initial  interaction)  then  the  maximum  number  of  nucleons 
which  could  be  in  statistical  equilibrium  when  this  ensemble 
emitted  an  a can  be  estimated  by  equating  2 (v^/v,^)  = (8/N*)  (q/P„) 
to  our  observed  a asymmetries,  where  N*  is  the  number  of  nucleons 
in  the  ensemble,  and  q is  the  momentum  transfer  of  the  virtual 
photon  which  initiated  the  reaction.  From  our  “‘Ag  and  *^Au  (e,a) 
mean  asymmetries  in  the  energy  interval  between  the  Coulomb  barrier 
and  25  MeV  (.16  and  .15,  respectively)  we  estimate  N*  < 18. 

The  intranuclear  cascade  model,  INC,  of  Bertini  ^ can 
roughly  account  for  both  the  inelastic  proton  energy  and  angular 
distributions  of  29  and  62  MeV  proton- induced  reactions^^.  Like  the 
pre-equilibrium  statistical  model,  the  basic  assumption  of  the  INC 
is  that  the  reaction  can  be  represented  by  a sequence  of  two-body 
interactions.  Unlike  the  pre-equilibrium  model,  the  INC  is  a 
microscopic  model  and  the  trajectory  of  a given  particle  is 
followed  through  the  nucleus  by  means  of  the  "correct"  free 
particle  scattering  cross  sections.  Bertini  ^ al^^  found  that  the 
fit  to  the  inelastic  (P/P')  data  of  Bertrand  and  Peele^^  could  be 
improved  if  their  calculation  included  four  particle  clusters  such 
as  a clusters. 

We  now  explore  a model  to  explain  our  (e,a)  data  in  which  the 
virtual  photon  is  absorbed  by  a quasideuteron.  In  this  model,  the 
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fast  (p,n)  pair  which  results  from  photodisintegration  of  the 
quasideuteron  then  interacts  with  the  remaining  nucleons  in  the 
nucleus  as  in  the  (p,a)  and  (n,Q:)  reactions.  This  model  is 
suggested  by  the  general  trends  of  the  (e,a)  cross  sections  shown 
in  Figs.  19  and  20  which  are  consistent  with  the  absorption  of  a 
virtual  photon  of  mean  energy  between  80  and  110  MeV.  It  is  also 
suggested  by  the  similarity  of  the  ratio  of  the  (p,a)  to  (p,p') 
inelastic  proton  scattering  cross  sections  integrated  over  the  same 
particle  energy  interval.  These  ratios  derived  from  Bertrand  and 
Peele's^^,  BP,  (p,a)  and  (p,p')  data  for  incident  proton  energies  of 
29  and  62  MeV  are  .1  and  .9  for  Ag,  and  .05  and  .09  for  Au. 
Since  quasideuteron  photon  absorption  at  Eq  = 110  MeV  would  result 
in  fast  neutrons  and  protons  with  energies  a between  20  and  50  MeV, 
we  should  compare  the  average  of  the  ratios  obtained  from  (p,a)  and 
(p,p')  at  29  and  62  MeV  of  .15  for  Ag  and  .07  for  Au  to  our  data. 
Our  ratios  are  .15  for  Ag  and  .045  for  Au  at  Eq  = 110  MeV. 


This  phenomenological  model  would  also  predict  that  the  (e,a) 
angular  distributions  would  be  approximated  by  folding  the  (e,p) 
and  (p,a)  angular  distributions.  Hence, 


I do 
\ dQ 


oe 


(cos0gp)  Pf,(cos0^J  dQp 


V'V'P, <COS0,J  / (2t  ^ 1) 


(18) 


where  A®’’,  and  A^",  are  the  coefficients  of  an  expansion  of  the  (e,p) 
and  (p,a)  angular  distributions  in  Legendre  polynomials.  Bertrand 
and  Peele's^^  (P/Cr)  and  our  (e,p)  angular  distribution  data  averaged 
over  20  < Tp  < 25  MeV,  results  in  dcj/dn)^^  a p^  + .45P,  + .OlPj  +.OIP3 
for  Tp  = 29  MeV  and  da/dn)^  « p^  + .32Pi  - .OOlPj  +.OIP3  for  Tp  = 62 
MeV.  Our  measured  da/dn)^^  a Pq  + .41Pi  - .24P2  +.OIP3.  If  we  use 
the  *^Au(p,a)  angular  distribution  data  of  Lefort  ^ al^^  at  a 
bombarding  proton  energy,  Tp,  of  157  MeV,  then  Eq.  (9)  predicts 
da/dfi)^  « nPo  + .14Pj  - .OlPj  +.OIP3,  which  does  not  agree  nearly  as 
well  with  our  ‘’^Au  da/dn)^  as  does  the  average  of  the  Tp  = 29  and 
62  MeV  data  of  BP. 

Lastly,  we  consider  the  possibility  that  a certain  fraction  of 
our  yields  might  arise  from  direct  interaction  between  an  a cluster 
and  a virtual  photon  in  the  sense  of  the  direct  photodisintegration 
models  developed  by  Courant^  and  others^^.  Such  models  have  been 
useful  in  understanding  portions  of  (7,p)  and  (7,n)  cross  sections 
and  angular  distributions  above  the  giant  resonance  in  light 
nuclei.  Since  we  only  detected  alphas  with  energies  < 26  MeV,  this 
model  might  be  expected  to  apply  to  the  asymmetric  (oc  sin^0cos0) 
component  of  our  da/dn  which  are  produced  by  the  absorption  of 
virtual  photons  with  energies  < 50  MeV.  At  higher  excitation 
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energies  a large  fraction  of  the  virtual  photon  energy  would  be 
retained  by  the  residual  nucleus,  which  is  not  compatible  with  the 
assumptions  of  this  model. 


If  a certain  fraction  of  our  asymmetric  alpha  and  proton 
yields  could  be  associated  with  a direct  photodisintegration 
process,  say  by  the  existence  of  a precise  correlation  between 
particle  kinetic  energy,  virtual  photon  excitation  energy,  and  the 
excitation  of  relatively  simple  configurations  in  the  residual 
nucleus,  these  yields  together  with  the  ratio  of  theoretical 
estimates  of  the  direct  photodisintegration  (7,0:)  and  (7,p)  cross 
sections  could  be  used  to  estimate  the  average  number  of  a clusters 
in  nuclei.  Carver^^  made  such  estimates  based  on  total  a yields 
from  V bombarded  by  32  MeV  bremsstrahlung,  assuming  only  El  photon 
absorption  for  both  alphas  and  protons.  Neither  our  data  nor 
previous  (7,  a)  data**’^^’^^  establishes  the  correlation  suggested  above 
and  therefore  the  principal  value  of  this  discussion  is  to 
establish  rough  ratios  of  A3/A0  for  the  (7,p)  and  (y,a)  angular 
distributions  which  would  be  expected  from  direct  photo- 
disintegration. In  spite  of  these  cautions,  we  will  obtain  an 
estimate  of  the  asymmetric  component  of  our  (e,p)  and  (e,a)  yields 
which  is  consistent  with  the  absorption  of  35  MeV  virtual  photons 
based  on  our  angular  distributions  and  excitation  data  for  Eq  = 50, 
75,  80,  110,  and  115  MeV.  We  then  assume  this  component  is 
produced  by  direct  photodisintegration  by  interference  between 
photons  of  El  and  E2  multipolarities.  Since  we  will  compare 
portions  of  the  a and  proton  energy  spectra  several  MeV  above  the 
a Coulomb  barrier,  we  assume  the  influence  of  the  barrier  in 
inhibiting  the  a yields  should  be  negligible.  If  (Z  - 2N„)  protons 
participate  in  the  direct  photodisintegration,  the  number  of  a 
clusters,  N„,  predicted  by  the  ratio  of  the  a to  proton  E1-E2 
interference  cross  sections  is 


« Z 


1 + 2iV/ 


where  iV,°  = 


(^qE1qE2)  ^ 


i£lO^ 


Atteiip 
At ten. ' 


( YS\ 


Asymmetric 


(19) 


and  Z is  the  nuclear  atomic  number.  Here,  e®'  and  e^  are  the 
effective  charges,  is  the  E1E2  interference  operator,  and 

and  Pp  are  the  density  of  states  for  the  a and  proton.  The  factor 
Ar  has  been  included  to  account  for  the  absorption  of  alphas  and 
protons  in  the  nucleus,  although  in  a strict  interpretation  of  the 
direct  photodisintegration  model,  this  factor  should  not  appear. 
We  use  our  ’Be  (e,a)  and  (e,p)  data  shown  in  Fig.  (14)  and  estimates 
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of  the  effective  number  of  a clusters  in  ^Be  obtained  from  (p,p'a) 
measurements  to  obtain  an  estimate  of  in  ^Be.  chant^^  estimates 
is  .451.15  from  the  analysis  of  ’Be(p,p'a)  data.  Based  on  our 
’Be  (e,a)  and  (e,p)  energy  spectra  at  34°  and  90°  obtained  with  Eq 
= 110  MeV,  Yr  for  alphas  and  protons  of  kinetic  energies  2014  MeV, 
corresponding  « to  the  absorption  of  3 8 MeV  photons  for  both  alphas 
and  protons,  is  .27.  Measurements  of  Buchnea  ^ show  the 

’Be(7,a)^He  reaction  is  ss  as  probable  as  the  ^Be  (7 , a)  n'’He  reaction. 
The  a+a+n  configuration  is  included  in  N^;  however,  since  this 
configuration  only  contributes  to  the  isotropic  yields  we  reduce 
our  estimate  of  by  . 5 . Using  Eg.  (19)  with  Ar  = and  = .22 
corresponding  to  = .06, 


(A-4)  ^ 


+ 8 (Z-2) 
+ iZ-l) 


~ .11. 


(20) 


This  value  for  M^  can  be  compared  with  a theoretical  estimate 
obtained  by  assuming  both  the  a and  proton  are  bound  in  a square 
well  in  an  relative  S state.  For  such  systems. 


M 


Be 


00 

[‘^pj 

(k„-S„)2cos62^ 

(21) 


where  S is  the  separation  energy,  and  P2„,  pjp  are  the  (p,*Li)  and 
(a,^He) , 1=2,  elastic  scattering  phase  shifts. 


To  evaluate  the  attenuation  factor,  Ar,  the  mean  free  paths, 
Xo  and  Xp  are  needed.  Igo  ^ al^  estimated  the  a mean  free  path  in 
nuclear  matter  of  5%  of  the  central  nucleon  density  to  be 

X,  ~ [ . 0215ln(r„)  + .318]  fm,  (22) 


from  data  for  T„  between  10  and  10^  MeV.  The  proton  mean  free  path 
obtained  from  an  optical  model  analysis  is^’ 


= 1.3 


1 + 


68.6 


fm, 


(23) 


We  estimate  the  attenuation  factors  for  a spherical  nucleus  with  a 
charge  density  parameterized  by  a Fermi  distribution'*®.  The  a 
distribution  was  assumed  to  be  a Gaussian  with  a maximum  located  at 
a distance  from  the  center  of  the  nucleus  suggested  by  Brink  and 
Castro*.  The  proton  distribution  was  taken  to  be  the  difference 
between  the  Fermi  and  a distributions.  The  relative  attenuation 
factor,  Ar  is  given  by 
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(Z-2JV;)F,(A,J 


(24) 


= 


This  result  is  derived  in  Appendix  II  and  values  of  F„(X„)  , F^(Xp)  , 
and  FpCXp)  are  given  in  Table  I.  For  a nucleus  of  mass  number  A, 


= 


8 (A  4 ) i[^e 
5(A-1)  ' 


(25) 


Values  of  = .08  for  T„  = 19  MeV  and  = .011  for  T„  = 23.5  MeV 
were  obtained  from  data  shown  in  Figs.  (15)  and  (16)  , and  values  of 
Ar  are  given  in  Table  II.  Finally  we  obtain  an  estimate  based  on 
the  direct  photodisintegration  model  of  N„  = 6±3  for  “'Ag  and  9±3  for 
'^Au. 


V.  CONCLUSIONS 

The  asymmetric  component  of  our  a-angular  distributions  cannot 
be  understood  in  terms  of  existing  statistical  models.  We  do  not 
know  of  completed  pre-equilibrium  or  precompound  model 
calculations'**  which  include  a rigorous  treatment  of  angular 
momentum.  While  such  calculations  are  underway'*^,  we  seriously 
doubt  that  calculations  which  treat  the  fabrication  of  an  alpha  by 
sequential  nucleon  pickup  can  simultaneously  reproduce  our  (e,p) 
and  (e,a)  angular  distributions.  Hence,  we  see  no  alternative  to 
attributing  the  asymmetric  component  of  our  yields  to  a direct  or 
semidirect  process.  A phenomenological  model  in  which  the  incident 
virtual  photon  is  absorbed  by  a quasideuteron^*  resulting  in  a fast 
(p,n)  pair  which  then  interacts  with  the  remaining  nucleons  via  the 
same  mechanisms  as  in  (p,a)  and  (n,a)  reactions  seems  to  best  fit 
our  data.  In  this  context,  the  arguments  used  to  relate  (p, a) 
and  (n,a)*®  reactions  to  nuclear  surface  clusters,  can  be  used  to 
relate  our  (e,a)  yields  to  nuclear  surface  a clusters. 
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TABLE  I.  Attenuation  factors  for  alphas  and  protons  in  Ag  and  Au.  The  sk: 
thickness  parameter  z of  the  Fermi  charge  distribution  = .527  fm.  Details  ( 
the  calculation  from  which  the  attenuation  factors  were  obtained  are  given  ; 
Appendix  II. 


c 

(fm) 

T 

(MeV) 

(fm) 

(fm) 

F.(X„) 

X 10’ 

F.(Xp) 

X 10 

Ff(>V 

X 

5 .20 

19 

6.78 

1.0 

. 37 

.64 

. 

5.20 

19 

6.78 

2 . 0 

. 37 

. 63 

. 

5 . 20 

19 

7 . 89 

2 . 0 

. 53 

.63 

. 

5.20 

19 

6.78 

. 1 

.23 

.64 

, 1 

5.20 

19 

7 . 89 

. 1 

1.10 

.64 

5.20 

30 

6.78 

1.0 

.40 

.44 

6.38 

23 . 5 

7 .96 

1.0 

.28 

.41 

6.38 

23.5 

7.96 

2 . 0 

.28 

.40 

6.38 

23 . 5 

8.96 

2 . 0 

.40 

.40 

6.38 

23 . 5 

7.96 

. 1 

. 18 

. 36 

6.38 

23.5 

8.96 

. 1 

.86 

.39 

6.38 

30 

7 . 96 

1.0 

.29 

.31 

i 


ABLE  II.  The  and  “‘Ag  (e,a)  and  (e,p)  angular  distribution  asymmetries 

^pressed  as  <asymmetry>  = ao  + a,<T>,  where  T is  the  kinetic  energy,  for  Eq 
50  and  110  MeV.  The  asymmetry  is  computed  over  the  range  of  T indicated, 
first  column  labeled  <asymmetry>)  and  over  an  energy  interval  10  MeV  wide 
nd  starting  at  the  T where  the  asymmetry  becomes  positive  (second  column 
abeled  <asymmetry>. 


Reaction 

Range 

(MeV) 

ao 

a, 

/MeV 

<asymmetry> 

Range 

<asymmetry> 

AT  = 10  MeV 

Eo 

(MeV) 

Au  ( e , a ) 

16<T„<26 

-.98±.21 

.061.01 

. 341.30 

. 321.29 

50 

Au  ( e , a ) 

16<T„<26 

-.771.08 

. 0441.004 

. 151.11 

.221.11 

110 

Au(e,p) 

8<Tp<2  3 

-.421.04 

. 051.01 

.361.16 

.251.14 

50 

Au(e,p) 

8<Tp<2  3 

-.351.02 

.0391.001 

.251.03 

.201.02 

110 

Ag(e,a) 

8<T„<2  6 

-.381.03 

.0321.002 

. 161.04 

. 161.05 

50 

Ag(e,a) 

8<T„<26 

-.551.04 

.0481.003 

.271.06 

.241.06 

110 

Ag(e,p) 

8<Tp<22 

-.221.02 

.031.01 

.231.15 

.151.13 

50 

8<T,<22 

-.261.05 

.04+. 01 

.341.16 

.201.13 

110 

l( 
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APPENDIX  I. 


Here  9?(q)  for  an  oscillator  well  {See  Eq.  (12)}  is  given. 

^(g)  = J , , where, 

^(g)  = Po^[llp/  + 4 (g+p^) 

- [ 1 1 + 4 {pj-  - p^g)  - 4p^^  ( 7p^2  + 4p^2 ) ] g -Pc^/Po^ 


9y/^gp^ 


erf 

' 2g+p^' 

- erf 

Pc  ] 

X ^Po  ; 

,y2Poj. 

(26) 


^(g)  = Po^{[llPo^+4  (g+p^.) 

- [llp/-Pe(g  + Pc)]e'^‘^'^‘^°'} 


+ 9sf2gpJ 


erf 


^ 2g+Pc'' 


( v^p 


+ erf 


^ Pc 


O / 


>/2P 


O/J 


The  quantities  p^,  q,  and  are  defined  in  the  text.  In  deriving 
9?(q)  , we  have  assumed  that  nucleon  momenta  p > Pc  can  be  neglected. 
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APPENDIX  II. 


We  calculated  the  a and  proton  attenuation  factors  by  assuming 
an  a particle  distribution, 


(27) 


and  normalized  so  that 

fp,(r}d^r  = 2N, 


(28) 


The  overall  charge  distribution 

Pf. 


(r-z) 

1 + e ^ 


Pf(^^  = 
fp,(r)d^r  = Z. 


was  given  by  a Fermi  distribution 
normalized  so  that 

(29) 


The  Fermi  distribution  parameters  c and  z were  talcen  from  Hahn  ^ 
aj..'*^  The  proton  charge  distribution  was  pp(r)-p„(r)  . The  parameter 

= [zJn(2)  + /2  (30) 


as  suggested  by  Brink  and  Castro^  and  a„  was  treated  as  a free 
parameter  subject  to  p“(r)  < p'’(r)  and  2N„  < Z.  The  a and  proton 
mean  free  paths  were  those  given  in  the  text.  The  attenuation 
factor,  A„,  for  alphas  was  then 


_lr-jd 

ff  e (r)  dMr) 

nuclear 

volume 

fj  p.(r)d^rdQ, 

nuclear 

volume 


(31) 


The  integration  is  over  the  volume  and  surface  of  a sphere  of 
radius  R such  that  Pp(R)  <<  .05pro* 


Ap(X^) 


ZF^(V  -2N^F^{X^) 
iZ-2Nj 


(32) 


The  dependence  of  the  attenuation  factors  on  r„  and  a„  is  explicitly 
displayed  in  Table  II. 
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FIGURE  CAPTIONS  I. 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Pulse  height  spectrum  produced  in  a focal  plane  detector 
by  the  electrodisintegration  of  ^Be  at  an  electron 
bombarding  energy  of  110  MeV.  The  kinetic  energy  of  the 
alphas  and  protons  is  4.5  MeV. 

‘^Au(e,a)  energy  distributions.  The  o is  for  an 
electron  bombarding  energy,  Eq  of  110  MeV.  The  □ is  for 
Eq  = 50  MeV.  The  energy  distributions  also  are  A^,. 

Ratio  of  coefficients  of  a least-squares  fit  of  ‘^Au(e,a) 
angular  distributions  for  Eq  = 110  MeV  to  an  expansion  in 
Legendre  polynomials  which  includes  terms  to  £ = 3 . 

Same  as  Fig.  3 except  that  Eq  = 50  MeV.  The  data  were  fit 
to  an  expansion  in  Legrendre  polynomials  which  included 
terms  to  £ = 2 . 


Fig.  5.  *’^Au(e,p)  energy  distributions.  The  o is  for  an  electron 

bombarding  energy,  Eq  of  110  MeV.  The  □ is  for  Eq  = 50 
MeV.  The  energy  distributions  also  are  A^. 

Fig.  6.  Ratio  of  coefficients  of  a least-squares  fit  of  ‘’^Au(e,p) 
angular  distributions  for  Eq  = 110  MeV  to  an  expansion  in 
Legendre  polynomials  which  includes  terms  to  £ = 3 . 


Fig.  7.  Same  as  Fig.  6 except  that  Eq  = 50  MeV.  The  data  were  fit 
to  an  expansion  in  Legendre  polynomials  which  included 
terms  to  £ = 2 . 


Fig.  8.  “‘Ag(e,a)  energy  distributions.  The  target  was  naturally 

occurring  Ag.  The  o is  for  Eq  = 110  MeV,  and  □ is  Eq  = 
50  MeV.  The  energy  distributions  also  represent  A^.  Also 
see  Fig.  3 for  ‘’^Au(e,a)  energy  distributions. 


Fig.  9.  Ratio  of  coefficients  of  a least-squares  fit  of  “‘Ag(e,a) 
angular  distributions  for  Eq  = 110  MeV  to  an  expansion  in 
Legendre  polynomials  which  includes  terms  to  £ = 3 . 

Fig.  10.  Same  as  Fig.  9 except  that  Eq  = 50  MeV.  The  data  were  fit 
to  an  expansion  in  Legendre  polynomials  which  included 
terms  to  £ = 2 . 


Fig.  11.  ""Ag(e,p)  energy  distributions.  The  target  was  naturally 

occurring  Ag.  The  o is  for  Eq  = 110  MeV,  and  □ is  for  Eq 
= 50  MeV.  The  energy  distributions  also  represent  Aq. 

Fig.  12.  Ratio  of  coefficients  of  a least-squares  fit  of  ^tAg(e,a) 
angular  distributions  for  Eq  = 110  MeV  to  an  expansion  in 
Legendre  polynomials  which  includes  terms  to  £ = 3 . 
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Fig.  13.  Same  as  Fig.  12  except  that  Eq  = 50  MeV.  The  data  were 
fit  to  an  expansion  in  Legendre  polynomials  which 
included  terms  to  £ = 2 . 


Fig.  14.  ’Be(e,a)  energy  distributions.  The  o is  for  Eq  = 110  MeV, 
and  the  □ is  for  Eq  = 50  MeV.  The  energy  distributions 
also  represent  Aq. 

Fig.  15.  *’^Au(e,a)  yields/MeV/sr/electron  as  a function  of  the 

bombarding  energy  Eq,  at  angles  between  the  direction  of 
the  emitted  a and  direction  of  the  incident  electron  beam 
of  34°  and  90°.  The  a are  for  an  a kinetic  energy,  T„  = 
16.0  ± .9  MeV,  o for  T„  = 17.9  ± .9  MeV,  <>  for  T„  = 19.8 
± .9  MeV,  V for  T„  = 21.7  ± .9  MeV,  □ for  T„  = 23.5  ± .9 
MeV. 


Fig.  16  “‘Ag(e,a)  yields/MeV/sr/electron  as  a function  of  the 

bombarding  energy  Eg,  at  angles  between  the  direction  of 
the  emitted  a and  direction  of  the  incident  electron  beam 
of  34°  and  90°.  The  a are  for  an  a kinetic  energy,  T„  = 
14.4  ± .9  MeV,  <>  for  T„  = 15.7  ± .9  MeV,  v for  T„  = 16.9 
± .9  MeV,  ■ for  T„  = 18.7  ± .9  MeV,  » for  T„  = 20.4  ± .9 
MeV. 


Fig.  17.  ^’^Au(e,p)  and  “*‘Ag(e,p)  yields/MeV/sr/electron  as  a 

function  of  the  electron  bombarding  energy.  Eg,  at  34°. 
The  upper  set  of  data  points  is  for  *’^Au(e,p)  and  here  o 
is  for  Tp  = 10.9  ± 1.4  MeV,  v is  for  T„  = 13.7  ± 1.4  MeV, 
A is  for  Tp  = 19.3  ± 1.4  MeV,  □ is  for  Tp  = 19.3  ± 1.4 
MeV,  <>  is  for  Tp  = 24.0  ± 1.4  MeV.  The  lower  set  of 
data  pints  is  for  “•*Ag(e,p)  and  here  o is  for  Tp  = 7.5  ± 
.9  MeV,  □ is  for  Tp  = 9.2  ± .9  MeV,  <>  is  for  Tp  = 10.9 
± .9  MeV,  • is  for  Tp  = 12.6  ± .9  MeV,  a is  for  Tp  = 18.2 
± 1.3  MeV  and  v if  for  Tp  =20.8  ± 1.3. 

Fig.  18.  ’Be(e,a)  yields/MeV/sr/electron  as  a function  of  the 
bombarding  electron  energy.  Eg,  at  34°.  The  o refer  to 
T„  = 6.7  ± 1.9  MeV,  □ refer  to  T„  = 8.5  ± 1.9  MeV,  v refer 
to  T„  = 10.4  ± 1.9  MeV,  <>  refer  to  T„  = 12.3  ± 1.9  MeV, 
A refer  to  T„  = 14.2  ± 1.9  MeV,  ★ refer  to  T„  = 16.0  ± 1.9 
MeV,  • refer  to  T„  = 14.2  ± 1.9  MeV. 


Fig.  19.  ‘’^Au  da(k)/dn  derived  from  excitation  data  shown  in  Figs. 

14  to  17.  The  curve  labeled  (a)  is  the  da(k)/dn  for 
protons  with  kinetic  energy  Tp  = 23  ± 1 MeV  at  0 = 34°. 
The  curves  labeled  (b)  and  (c)  are  da(k)/dn  for  alphas  at 
90°  and  34°,  respectively. 


Fig.  20.  “‘Ag  (e,p)  and  (e,cr)  cross  sections.  The  □ is  the  proton 
da(k)/dn,  and  o is  the  alpha  da(k)/dn  at  the  indicated  a 
kinetic  energies  for  6^  =90°. 
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Fig.  21.  The  number  of  recoil  alphas  with  energies  above  a cutoff 
kinetic  energy,  T^,  from  elastic  '*He(e,e')a  as  a function 
of  the  incident  electron  energy,  Eq. 


Fig.  22.  The  number  of  recoil  prptons  with  energies  above  a cutoff 
kinetic  energy,  T^,  from  elastic  '‘H(e,e')p  as  a function 
of  the  incident  electron  energy,  Eq. 

Fig.  23.  Comparison  of  ^^Au(e,p)  energy  distributions  for  Eq  =50 
(□)  and  110  MeV  (O)  with  the  predictions  of  the 
statistical  model  described  in  the  text.  The  theoretical 
prediction  for  the  peak  of  the  energy  distribution  had  to 
multiplied  by  4.26  to  bring  the  prediction  into  agreement 
with  experiment. 


Fig.  24.  Comparison  of  ‘’^Au(e,a)  energy  spectra  for  Eq  =50  (□)  and 
110  MeV  (O)  with  the  predictions  of  the  statistical  model 
described  in  the  text  and  with  the  same  input  data  and 
parameters  used  in  the  calculation  for  the  ^^Au(e,p) 
cross  section  shown  in  Fig.  23.  In  this  case,  the 
theoretical  prediction  had  to  be  multiplied  by  3.1  x 10*^ 
to  bring  the  peak  prediction  into  agreement. 


Fig.  25.  Comparison  of  a/eq.  quantum  obtained  from  this  work  o and 
from  500  MeV  bremsstrahlung,  histogram.  Ref.  16.  The 
smooth  curve  is  a calculation  taken  from  Ref.  16  and 
described  in  the  text. 
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